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ABSTRACT 

The Herschel Space Observatory has had a tremendous impact on the study of 
extragalactic dust. Specifically, early-type galaxies (ETG) have been the focus of sev¬ 
eral studies. In this paper we combine results from two Herschel studies - a Virgo 
cluster study HeViCS and a broader, low-redshift H-ATLAS/GAMA study - and con¬ 
trast the dust and associated properties for similar mass galaxies. This comparison 
is motivated by differences in results exhibited between multiple Herschel studies of 
early-type galaxies. A comparison between consistent modified blackbody derived dust 
mass is carried out, revealing strong differences between the two samples in both dust 
mass and dust-to-stellar mass ratio. In particular, the HeViCS sample lacks massive 
ETG with as high a specific dust content as found in H-ATLAS. This is most likely 
connected with the difference in environment for the two samples. We calculate near¬ 
est neighbour environment densities in a consistent way, showing that H-ATLAS ETG 
occupy sparser regions of the local Universe, whereas HeViCS ETG occupy dense re¬ 
gions. This is also true for ETG that are not iJersc/ieZ-detected but are in the Virgo 
and GAMA parent samples. Spectral energy distributions are fit to the panchromatic 
data. From these we find that in H-ATLAS the specific star formation rate anticor¬ 
relates with stellar mass and reaches values as high as in our Galaxy. On the other 
hand HeViCS ETG appear to have little star formation. Based on the trends found 
here, H-ATLAS ETG are thought to have more extended star formation histories and 
a younger stellar population than HeViCS ETG. 

Key words: methods: statistical - astronomical data bases: surveys - galaxies: ellip¬ 
tical and lenticular, cD - galaxies: evolution - submillimetre: galaxies 
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1 INTRODUCTION 


Dust is a fundamental component of the interstellar medium 
(ISM) of galaxies for the thermodynamics and chemistry of 
the gas, for the dynamics of the accretion in dense star¬ 
forming clouds, and for the attenuation of UV/blu e radia¬ 
tion and its re-emission in the far infra-red (FIR; IPraind 
l2003h . The relative amount of dust varies strongly with 
galaxy type, increasing by about three orders of magnitude 
on ave rage along the Hubble sequence (e.g. ICortese et alJ 
l2012bl l. The connection between dust and chemical evolu¬ 
tion also varies with galaxy type. In late-type galaxies (LTG) 
dust is strongly linked with star formation (SF), both be¬ 
cause it serves as a catalyst for the formation of molecular 
gas necessary for SF and because, being heated mostly by 
young stars, its emission traces the regions of SF. The same 
paradigm does not necessarily apply to early-type galaxies 
(ETG; comprising of ellipticals and lenticulars), where dust 
can be heated by the radiation field produced by evolved 
stars and can be more diffuse, therefore not serving as a SF 
catalyst. In addition, ETG, particularly those in clusters, 
can have much larger amounts of hot gas than LTG, not 
favouring the presence of dust. 


It is therefore important to study dust in ETG 
separate ly from LTG and the Herschel Space Obser¬ 
vatory { Pilbratt et ^ I I2OIOII has al lowed several such 
detai l ed studies iRowlands et al. |2012| : ISmith et al.l 
2 OI 2 I: di SeregQ Alighieri et al.l l2013l. S13 hereafter, and 


Agius et al.ll2013l . A13 hereafter). In particular S13 searched 


for dust in a large optical sample of 910 ETG in the Virgo 
cluster, extendi n g als o to dwarf ET G, using the PACS 
JPoglitsch et al.l I2OI0II and SPIRE (iGrifEn et al.l I2OI0I I 
instruments, and found it in 17% of the elliptical galaxies, 
in 40% of lenticulars (SO -I- SOa) and in about 3% of the 
dwarfs (dE -|- dSO). They showed that the presence of dust 
does not correlate with the presence of neutral gas (HI) and 
the dusty ETG do not appear to have bluer B-H colours, 
i.e. to be more star-forming than the non-dusty ones. On 
the other hand A13 searched for dust, also with PACS and 
SPIRE, in a sample of 771 brighter ETG (Mr < —17.4mag) 
over a very large volume (144 deg^ and 0.013<z<0.06) 
and found it in 20% of the ellipticals and in 38% of the 
lenticulars. Not only are these detection rates high, but 
also the relative amount of dust is higher than in HeViCS 
and the dusty ETG have bluer colours, suggesting that 
they may be forming more stars. Furthermore S13 found 
a dependence of the dust temperature on the stellar mass 
and on the average B-band surface brightness within the 
effective radius, but A13 did not. 


These differences may be an effect of the environment, 
since H-ATLAS covers a wide range of environments, po¬ 
tentially hosting extreme mergers or interactions, whilst 
HeViCS is limited to a clu ster environm e nt, w hich may 
be less favourable for dust. iBourne et al.l ll2012li pointed 
out that environment was a possible influence on the rel¬ 
atively low dust masses in Virgo cluster ETG, in compari¬ 
son with their results for stacked red sequence galaxies from 
the GAMA survey. Alternatively the dusty H-ATLAS ETG 
could represent younger versions of the ‘standard’ ETG; 
i.e. they may have formed recently, or have more extended 
star formation histories. In fact, if galaxy evolution is in¬ 
fluenced by the environment, ETG in sparse environments 


are more likely to be at an earlier stage i n their evolu¬ 
tion than ETG in dense environments (e.g. iThomas et al.l 
I 2 OO 5 II . Some differences may also be explained by the dif¬ 
ferent models used to estimate the dust mass or by limita¬ 
tions in the models or wavebands used. Although a single 
modihed blackbody (modBB) fitting approach gives a good 
estimate of the mass of cold ^ diffuse dust grain s in the ISM 
of ga l axies at all redsh ifts iD_unng_et_alJ 2000|;_ Blain_et al.l 
I 2 OO 3 I : iPope et ^ I 2 OO 6 I : iDve et al.l l20ld : Bianchil |2013|) . it 


may not account for the emission from dust in warmer me¬ 
dia, such as the grains surrounding birth clouds of hot, 
young stars, or colder dust that needs multiple tempera¬ 
ture models. The addition of further modBBs allowing for 
diffe rent temperatures have been shown to improve such 
hts I Dunne fc Baled 1200 ll : ICalametz et aLll^lluDale et al.l 
I2OI2 1. Furthermore, given the wealth of panchromatic data 
for HeViCS and H-ATLAS/GAMA ETG, it is possible to ex¬ 
ploit multi-wavelength SED fits which consider stellar emis¬ 
sion at UV/optical wavelengths, the attenuation by dust and 
resultant emission in the infrared. 

Temi et al. (2009a,b) and Amblard et al. (2013) have 
further investigated the diversity of ETG by studying the 
physical properties of a sample of local E and SO galaxies. 
They find that many local SO galaxies are quite distinct from 
ellipticals, containing dust and cold gas in amounts that 
may be sufficient to generate appreciable star formation at 
rates as large as several M0/yr. However in this paper we 
cannot investigate in detail the differences between E and 
SO galaxies, since they are difficult to distinguish in the H- 
ATLAS/GAMA sample because of their distance and the 
limited spatial resolution of SDSS images. 

In this paper we compare the properties of the dust 
in ETG from both the HeViCS sample of S13 and the H- 
ATLAS/GAMA sample of A13, by characterising them in a 
uniform way and taking into account data at shorter wave¬ 
lengths, in an attempt to understand the differences and the 
reasons causing them. 

This paper is laid out in the following manner. In Sec¬ 
tion [2] we describe the data available for the HeViCS and 
the H-ATLAS/GAMA ETG samples. In Section [3] we com¬ 
pare the two samples, while in Section |4] we discuss effects 
of environment. Panchromatic fits to the ETG SED are de¬ 
scribed in Section [S] and a discussion of derived parameters 
is given in Section (6] and is followed by our conclusions in 
Section 7. The photometric data obtained from WISE are 
described in Appendix A. We assume a flat Universe with 
Dm = 0.3, Ha = 0.7 and Ho = 70 km/s/Mpc. 


2 DATA 

This section summarises the different surveys utilised in con¬ 
structing the ETG samples that will be compared in this 
work. 

2.1 Far-Infrared Surveys 

The Herschel-Astrophysical Terahertz Large Area Survey 
(H-ATLA^l) is the widest open-time extragalactic survey 


^ http://www.h-atlas.org 
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with the Herschel Space Observatory, with one of its ulti¬ 
mate science aims being the investigatio n of the dust con¬ 
tent of the nearb y Universe at z<0.5 (jEales et al.l l201fll : 
iDunne et al.l I2OIII I . This survey samples over ~570 deg^ 
of sky, covering a range of environments in a uniform 
way. Their data collection process involved parallel imag¬ 
ing wit h Herschel’s two pho tometers, PACS (100/im and 
160^m; IPoglitsch et al. 2010ll a nd SPIRE (250^m, 350/im 
and 500^m: loriflin et al.ll201ol i. with a 60”s“^ scan rate. 
H-ATL AS has a 5g sensi tivity limit of 33.5 mjy/beam at 
250/rm dRigbv et al.ll201ll l - this corresponds to a dust mass 
of for a range of temperatures (15-30 K) at low 

redshift (z<0.06). 

H-ATLA S catalogues were co nstru cted from maps a s 
described in IPascale et al.l (120111 1 and llbar et ^ l|2010h . 
Source extraction was based on emission greater than 5 (t 
in any of the 3 SPIRE wavebands, d escribed in d e tail fo r 
the Science Demo nstration Phase in iRigbv et al.l ll201lll . 
ISmith et al.l l|201ll l gives a description of the likelihood-ratio 
analysis performed to identify r-band optical counterparts 
to the SPIRE sub-mm selected sources. Based on the re¬ 
sultant positional and photometric information for the in¬ 
dividual sources, PACS flux densities are measured using 
circular apertures placed at the SPIRE positions. Details of 
the Phase 1 H-ATLAS and the GAMA catalogues used in 
this paper can be found in A13. 

The Hersche l Virgo Clu s ter Survey (HeViCS 
iDavies et al. ] I 2 OIOI: iDavies et al.1 l2012ll is an audit of a 
large fraction (84 square degrees) of the Virgo Cluster 
in the same five Herschel bands as the H-ATLAS survey. 
This specifically samples the dense environment of a 
nearby cluster, going down to fainter luminosities than 
H-ATLAS. Additionally their observations are deeper than 
the H-ATLAS observations, with four linked cross-scans 
for HeViCS compared to a single cross-scan for H-ATLAS. 
HeViCS observations were performed in fast-parallel mode 
with PACS and SPIRE, with a scan rate of 60"s“^. The 
HeViCS 5cr sensitivity at 250Atm is 25-33 mjy for sources 
smaller than the PSF (lAuld et al.l 1201^ : S13); depending 
on the dust temperature this corresponds to a dust mass of 
<^0.2-1 X 10®Mq at the 17 Mpc distance of the main Virgo 
cluster cloud. A detailed account of the data collection, re¬ 


self-consis tent using an a perture matching technique de¬ 
scribed in iHill et al.l (l201llL 

The spectroscopic element of GAMA has just been 
completed at the Anglo-Australian Telescope (AAT), with 
the most recent Data Release (DR2) from GAMA giv- 
ing access to 70,000 new redshifts in the GAMA I regions 
d Hopkins et al.l l2013l : iLiske et al.l lin preoh . These are three 
regions of 48 square degrees each, centred at 9, 12 and 14.5 
hours (G09, G12 and G15) on the celestial equator. Spectro¬ 
scopic completeness limits are given as apparent petrosian 
magnitude ip^t <19.4 mag in the G09 and G15 fields, and 
rpet <19.8 mag in G12 for GAMA 1. 

The H-ATLAS equatorial fields coincide with those of 
GAMA, and matching between the two sets of data revealed 
'^10,000 counterparts, using the likelihood ratio method 
(Bourne et al. in prep). It is from within these counterparts 
that the H-ATLAS ETG sample is constructed, as described 
in A13 and H3.ll Section 2 of A13 gives details of the GAMAI 
databases used in this current paper. 

A very large set of data is available for galaxies in the 
Virgo cluster. In this nearby cluster we can obtain accu¬ 
rate galaxy morphological classifications and observe galax¬ 
ies covering a wide range of luminosities. The main orig¬ 
inal s ource of information is th e Virgo Gluster Catalogue 
('VCC. lBinggeli et al.ll98^.l l993l l which, together with Virgo 
SDSS data ( Davies et al. 2014l L remains the most complete 
optical catalogue, until the catalogu e of the Next Gener - 
ation Virgo Gluster Survey fNGVS. iFerrarese et al.ll2012h 
becomes publicly available. The VGC is complete to pho¬ 
tographic magnitude mpg = 18.0, but also contains fainter 
galaxies. GOLDMine (iGavazzi et al.l[2003l l provides a com¬ 
pilation of data on VCG galaxies. Useful additions are the 
GALEX Ultraviole t Virgo Gluster Survey papers (GUViCS, 
|Boselh__efL_alJ|201l|;LBi3selIj|| ^lj) a nd the HI survey for ETG 
of Idi SeregQ Alighieri et ah (20^. 


duct io n and flux measurements can be found in lAuld et al.l 
(l2013li and in S13. 


3 OVERVIEW OF ETG SAMPLES 

This section compares and contrasts two ETG samples, 
which are described below. Particular emphasis is placed on 
the differences between the classification criteria for these 
samples. A summary of the results from their parent papers 
is also given. 


2.2 Multi-Waveband Data 


3.1 H-ATLAS Sample 


The Galaxy and Mass Assembly (GAMAQ) is a spectro¬ 
scopic and photometric survey dedicated to constructing 
a galactic database which spans the electrom agnetic spec¬ 
trum from ultraviolet to radio wavebands ([Driver et al.l 
[ioil). This campaign is being supplemented by imaging 
from survey s such as the Sloan Digital Sky Survey in 
the optical (lAb azaiian et al.l l2 009ll . GALEX in the UV 
( Bianchi fc GALEX TeamI 1 19991) , UKI DSS-LAS in the N IR 
( Lawrence et al.ll2007l l and H-ATLAS (lEales et al.ll201ol ’) in 
the FIR/sub-mm; all these surveys have overlapping data 
within the same regions and the photometry has been made 


^ http://wiki.arcetri.astro.it/bin/view/HeViCS/WebHome 
® http://www.gama-survey.org 


From the sample of 771 ETG with M^ <-17.4 in the GAMA 
equatorial fields, the H-ATLAS detections comprise 220 Es 
and SOs with 5cr 250/im from Herschel. Thus there are 551 
ETG in our GAMA sample that are below this detection 
limit (the undetected H-ATLAS ETG). The optical coun¬ 
terparts in GAMA have reliability of association > 0.8. The 
morphological classification process for all these ETG is fully 
described in A13; briefly, it was based on visual classifica¬ 
tion of blue, green and red optical galaxy cutouts into six 
groupings of E, SO , SBOa, Sbc, SBbc and Sd galaxies (see 
I Kelvin et al.l f20l3 for a full account of this process). The 
galaxies classified in this way are GAMA I galaxies within a 
redshift range of 0.013<z<0.06 and complete to an absolute 
magnitude cutoff of M^ <-17.4 - these limits are therefore 
also applicable to the ETG sample. We note that this sample 
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lacks the faintest ETG dose to the Mr=-17.4 limit, proba¬ 
bly because these are dominated by the PSF and therefore 
excluded since they cannot be reliably classified (see above). 
We will discuss the effects of this selection in §[331 

The ETG sample was constructed from the H-ATLAS 
detected E and SO (which include both SO, SOa and Sa galax¬ 
ies) galaxies from within this classified set of galaxies, with 
additional criteria imposed to remove any potential spiral 
structure, edge-on disks, or small objects which may be 
barely-resolved and thereby possibly misclassified. Galax¬ 
ies with AGN and LINER si gnatures in the optical BPT 
diagram ([Baldwin et aljriOSli l were also removed, so as to 
only consider galaxies with a FIR/sub-mm SED dominated 
by thermal dust emission. This resultant sample contains 
73 Es and 147 SOs, a few examples of which are shown as 
g-band SDSS cutouts in Fig. [T] Please see A13 for details 
of comparisons of properties between Hubble types, for the 
H-ATLAS sample. 

A13 showed that the H-ATLAS sub-mm detected ETG 
sample had unusual characteristics in comparison to the un¬ 
detected ETG. In particular, both optical and UV-optical 
colours were typically quite blue, indicatin g some ongoing 
star formation in these systems (see also [Rowlands et al.l 
I2ni2h . The galaxy light profiles showed more exponential (or 
less centrally concentrated) luminosity distributions, which 
might indicate some recent merging activity, or may be an 
effect of dust attenuation, or a faint disk. Finally, an investi¬ 
gation of nearest neighbour galaxy surface density revealed 
that these ETG inhabit sparser environments than the non- 
detected ETG. 

ModBBs with emissivity spectral index (5= 2 and 350/xm 
mass absorption coefficient K35 o=4.54 cm^g“^ ([Punne et al.l 
l201lh were fit to the PAGS (100 and 160/im) and SPIRE 
(250, 350 and 500^m) data for 188 galaxies (see A13). 
These parameters are fixed to these values throughout this 
paper. A13 report a range of rest-frame dust tempera¬ 
ture between 9-30 K and a range of dust masses between 
8.1 X 10^-3.5x 10 ®Mq, with a mean dust-to-stellar mass ratio 
of logio(Md/M,)=-3.37. These results are the key parame¬ 
ters which will be investigated within this paper, in compar¬ 
ison to the nearby ETG in the Virgo Cluster. A summary of 
the characteristics of H-ATLAS and HeViCS samples (de¬ 
scribed in 113.21) is given in Table [T] 


3.2 HeViCS Sample 

For HeViCS, S13 utilise an input optical sample from the 
VCC, constr ained by ETG mor phology (as compiled in 
GOLDMine; ICavazzi et al.l l2003f l but not limited in any 
other respect. Therefore their sample of 910 ETG spans a 
range of magnitudes, from dwarf systems to the most mas¬ 
sive ETG, and contains classifications equal or earlier than 
SOa-SO/Sa types. 

S13 found 52 ETG by searching for Herschel counter¬ 
parts within one pixel (6") and signal-to-noise (S/N) greater 
than 5, in the parent sample described above. The reliabil¬ 
ity of these counterparts is fully discussed in S13. In order 
to make a fair comparison with the brighter galaxies in the 
GAMA/H-ATLAS sample, we selected a bright subsample 
of the 910 ETG of the input HeViCS sample, by applying a 
cutoff of SDSS Mr <-17.4, as with the GAMA galaxies. This 
subsample is complete for the Virgo Cluster and contains 


Parameter 

H-ATLAS Det. 

HeViCS Det. 


min 

max 

min 

max 

Sample Size (gals) 


220 


33 

Distance^, (Mpc) 

57.2 

265.4 

17.0 

32.0 

logic (M*) (Mq) 

8.9 

11.4 

8.7 

11.4 

Mr (mag) 

-18.2 

-23.1 

-17.4 

-23.1 

mr (mag) 

17.7 

13.3 

14.1 

8.1 

i^NUV (mag) 

22.7 

16.6 

18.1 

13.8 

F 250 (Jy) 

0.033 

0.770 

0.013 

7.992 

L 250 (W Hz-1) 

1 .6x10^2 

4.3x102-1 

4.4x102° 

2.8x1023 

logic(Md) (M 0 ) 

5.91 

8.54 

4.48 

6.67 

logic(Md/M,) 

-4.44 

-2.13 

-6.29 

-3.07 

Sgal (gals Mpc-2) 

0.001 

37.08 

29.19 

463.10 


Table 1. Parameters indicating the types of ETG found in the 
sub-mm detected H-ATLAS and HeViCS samples. Note that pa¬ 
rameters for the HeViCS sample only include the 33 ETG with 
Mr <-17.4. The parameters shown include sample sizes, ranges of 
luminosity distances, stellar masses, r-band absolute magnitude, 
r and NUV band apparent magnitudes, total 250 /rm fluxes and 
luminosities. Dust mass, dust-to-stellar mass ratio and environ¬ 
mental density ranges are also shown, calculated as described in 
the main text. 


123 ETG. Out of these, 33 ETG are detected at 250 /rm: 
their properties are given and contrasted to the H-ATLAS 
sample in Table [T] Sixteen of the nineteen ETG detected at 
250 /im, which have been removed, have formal GOLDMine 
dwarf classihcations. From here on, these 33 ETG will form 
the HeViCS detected sample. Correspondingly 90 Virgo 
ETG with Mr <-17.4 are undetected at 250/tm. The mag¬ 
nitude cutoff Mr <-17.4 also has the effect of removing all 
those ETG identified as possible contaminating background 
sources in S13. Therefore all the 33 HeViCS ETG considered 
here have secure identifications. Fig. [2]shows a few examples 
of Virgo ETG detected in HeViCS. VCC763 and VCC1535 
are shown to have centrally distributed sub-mm emission, 
but VCC881 is a special case, as the sub-mm contours are 
quite offset from the galactic centre. This is likely due to 
the streams linked wi th this galaxy, indicative of dwarf com¬ 
panion stripping (e . g. pTrin chieri fc di S ereeo Ahghieri|[l99l[ : 
ikennev et al.l[20^ : ljanowiecki et al.ll201C L 

Apparent r-band magnitudes fo r this ETG sample hav e 
been obtained from the work by [Cortese et al.l (l2012a[ L 
where they calculated UV and optical asymptotic mag¬ 
nitudes for the HRS galaxies, some of which are in the 
Virgo cluster. This provided AB r-band magnitudes for 148 
HeViCS galaxies; the remaining HeViCS galaxies have mag¬ 
nitudes calculated from the combination of B-band magni¬ 
tudes from GOLDMine and the average (B-r) = 1.02±0.26 
colour obtained from these 148 galaxies. These apparent 
magnitudes are then converted to absolute magnitudes us¬ 
ing the GOLDMine distances and the appropriate galactic 
absorption. 

Stellar mas ses were estimated for these galaxies using 
the method of IZibetti et al.l (l2009ll whereby optical (and 
NIR when available) photometry and synthetic libraries are 
compared. Dust temperatures and masses were derived from 
modBB fits to the FIR/sub-mm data. For the 33 massive 
(i.e. non-dwarf) ETG and the same values for g and ftsso 
as the H-ATLAS sample, dust temperatures and masses are 
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Figure 1. Example images of sub-mm detected H-ATLAS galaxies with E (top) and SO (bottom) classifications. The images are 40” by 
40” SDSS g-band images with superimposed H-ATLAS 250/rm contours in red. These contour levels represent ~15, 35, 55, 80 and 100% 
of the 250/rm flux. Galaxy classification, catalogue ID and absolute r-band magnitudes are shown on the top-left of the images. The 18” 
SPIRE FWHM PSF, at 250/rm, is also shown in these images. 




Figure 2. Example images of sub-mm detected HeViCS galaxies with a variety of classifications. The images are 6’ by 6’ SDSS g-band 
images with superimposed HeViCS 250gm contours in red. These contour levels represent the following percentages of the 250gm flux: 
14, 43, 71 and 100% (VCC763), 18, 36, 55, 73 and 100% (VCC881), 3, 28, 52, 76 and 100% (VCC1535). Galaxy GOLDMine classification, 
identification and absolute r-band magnitudes are shown on the top-left of the images. Note that VGC 763 has a synchrotron component. 
The 18" SPIRE FWHM PSF is also shown in these images. 


given as 14.6-30.9 K and 3.0x 10‘^-4.7x IO^Mq. Their mean 3.3 Sample Comparison 
dust-to-stellar mass ratio is logio(Md/M*)=-3.93. 

There are some clear differences between the two samples 
which need to be addressed before proceeding with a com¬ 
parison of their properties. A primary concern is the dif¬ 
ference in galaxy distance (see Table [1} : HeViCS ETG are 
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located in the nearby Universe at a distance between 17 and 
32 Mpc, whereas H-ATLAS ETG are further away within 
a distance range of 56< (11 <269 Mpc with an average dis¬ 
tance of ~195 Mpc. As a consequence, H-ATLAS ETG will 
have lower spatial resolution, larger luminosity at the op¬ 
tical detection threshold, and a higher dust-mass detection 
threshold. HeViCS ETG are very well resolved and have 
lower detection thresholds at all wavebands. For H-ATLAS, 
this results in the morphological classihcation not being as 
detailed as that completed for HeViCS. Therefore H-ATLAS 
ETG can be identified as either E or SO galaxies, but we can¬ 
not distinguish safely between these two classes, nor detect 
any dwarf galaxies, which in any case are excluded by the 
Mr <-17.4 limit. Also the completeness at this limit of H- 
ATLAS ETG is not as good as for the HeViCS ETG for the 
reasons explained in Section 3.1. 

Given that this work will contain a statistical analysis of 
the properties of the two samples, it is important to consider 
whether the ETG sample sizes are statistically significant. 
Additionally, when comparing properties for the samples us¬ 
ing statistical tests, it is preferable for the sample sizes to 
be of similar orders of magnitude in order to obtain a fair 
analysis. The H-ATLAS sample contains 220 ETG detected 
in the FIR whereas HeViCS contains 33 ETG (within the 
H-ATLAS magnitude cutoff of Mr <-17.4). Both samples 
are large enough to run a Kolmogorov-Smirnov (KS) test, 
which is sensitive to fairly small differences between modest 
sized samples, to check whether the populations are similar. 
The sample sizes themselves are different with the HeViCS 
sample only containing ~15% of the H-ATLAS numbers, 
but the differences are not so large that such a test would 
be invalid. 

HeViCS dust masses and specihc dust masses are lower 
than those of the H-ATLAS sample by approximately an 
order of magnitude, even though the morphologies of the 
galaxies are similar and the stellar mass ranges overlap. This 
is due to the fact that the closest H-ATLAS ETG (those 
at Z'~0.013 or distance Al ~56 Mpc) are more than three 
times further away from us than the main Virgo cloud at 
17 Mpc, and on average the H-ATLAS ETG are still much 
further away (see Table [T]). Therefore smaller quantities of 
dust (at least ten times smaller) can be detected in HeViCS 
ETG than in H-ATLAS. These differences need to be taken 
into account to avoid possible biased conclusions about the 
properties of ETG as a class. 

Dust appears to be much more concentrated than stars 
in Virgo ETG , and more lumino us ETG have higher dust 
temperatures dSmith et al.ll20f^ .S13L The dust mass does 
not correlate clearly with stellar mass, while the dust-to- 
stellar mass ratio anticorrelates with galaxy luminosity. The 
dusty ETG appear to prefer the densest regions of the Virgo 
cluster. Contrary to H-ATLAS/GAMA ETG, the HeViCS 
ETG detected at 250jLtm are not bluer t han the undetected 
HeViCS ETG ddi Seregp Alighierill2013r) . 

The significant difference in 250^m luminosities be¬ 
tween the two samples also needs to be considered, es¬ 
pecially given the similarity in optical luminosity between 
the two samples. We calculate this parameter for all the 
ETG. H-ATLAS ETG have a range of luminosities from 
1.6x 10^^-4.3x 10^^ W Hz“^, whereas HeViCS have a range 
from 4.4xl0^°-2.76x10^® W Hz“^. We define a threshold 
luminosity, defined by the H-ATLAS flux limit of 33.5 mjy 


and maximum distance (z=0.06 ^ di:,=269 Mpc), which 
equates to 2.77x10^® W Hz“^. There are 52 ETG (24%) 
within H-ATLAS with detections below this threshold lu¬ 
minosity; therefore there is some overlap in FIR luminosity 
between the two samples. It is possible that the difference 
in luminosities is a direct result of the much larger volume 
covered by H-ATLAS/GAMA, which gives a larger chance 
of seeing rarer objects. 


4 EXPLORING ENVIRONMENTS 

Galaxies in the HeViCS sample are, by design, located in 
the Virgo Cluster. Although the cluster is an overall dense 
environment, the density is not homogeneous and will vary 
with position throughout the cluster. Conversely, because 
the ETG in the H-ATLAS sample were taken from a wide 
area of sky over a volume of redshift, they most likely belong 
to a range of environments, thus reaching lower densities. 
Therefore, a quantitative comparison of the environments 
inhabited by the ETG in both samples requires some form 
of estimation of the environmental density performed in a 
consistent manner. 

4.1 Nearest Neighbour Densities 

To calculate environmental densities we utilise nearest 
neighbour surface densities. This has a lready been done to 
some extent for the H-ATLAS sample dBroueh et al.llioisl . 
A13), although a bright magnitude limit of Mr <-20 was im¬ 
posed that may not accurately sample the true densities of 
these ETG. Nearest neighbour densities are now calculated 
which do not incorporate so bright a magnitude limit. 

Chris Beaumont’s IDL librarj0 is used to calculate a 
smoothed map of the coordinates of all the galaxies in the 
HeViCS and H- ATLAS sample regions r espectively, based 
on the method in iGutermuth et al.l (l2005l l. Although we are 
only interested in the environments of the sub-mm detected 
ETG, it is necessary to perform this routine on the entire 
population of galaxies within these regions to accurately de¬ 
pict the true density; this is the density-dehning population 
(DDP). For every galaxy, the algorithm calculates the dis¬ 
tance Djv to the N-th closest object and thus the surface 
density 

E,. = (1) 

The value of N chosen for these calculations is hve as this 
is a good compromise given the effects of survey edges. Ad¬ 
ditional calculations are then required to convert the units 
of the surface densities from objects per square degrees to 
objects per square Mpc. For the HeViCS galaxies, this is a 
straightforward conversion using the distance of the galaxy. 

These calculations are not as simple for the H-ATLAS 
sample; because of the large redshift range of the galaxies, 
iBrough et al.l (l2013l l and A13 limited the DDP for each sam¬ 
ple galaxy to a velocity cylinder ±1000 km s“^ over which 
the surface density is calculated, so that the latter is not in¬ 
fluenced by galaxies at large distances, which clearly cannot 

http://www.ifa.hawaii.edu/users/beaumont/code/ 
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Egai (Galaxies Mpc") - This Work 

Figure 3. A comparison of the nearest neighbour densities for the 
H-ATLAS sample as calculated in this work with those calculated 
in A13. The key difference in the calculation of these densities 
is the magnitude limit of the DDPs required to calculate these 
values. A one-to-one correlation is shown to aid comparison. The 
inset plot shows the quantitative difference between the two in 
histogram form. 


have any environmental effect. This is also repeated here. 
Once each DDP has been created, the procednre described 
above for the HeViCS densities can be rnn, and surface den¬ 
sities created. 

An additional restriction for these calculations is the 
imposed magnitude limit on the galaxies used to create the 
DDP. The Virgo Cluster galaxies can be detected down 
to much fainter magnitude limits than the higher redshift 
galaxies, and therefore surface densities for the latter are 
likely to be underestimated because dwarf galaxies which 
are detected in the Virgo cluster cannot be detected at 
higher redshifts. To avoid this, we set a magnitude limit 
on the DDP of Mr <-17.4 for both samples. This is the 
faintest limit which still ensures completeness also on the 
H-ATLAS/GAMA sample at 0.013 < a < 0.06. 

We test the calculated s urface densities by c omparing 
them with those derived by iBroueh et al.l (l2013l l and dis¬ 
played in A13. Fig. [3] shows a direct comparison between 
the two parameters, with a mostly linear relation defined. 
There are some galaxies from A13 which have higher surface 
densities than those calculated here: these are upper limits 
calculated for these galaxies flagged as lying at the edge of 
the GAMA regions. More importantly, the surface densities 
calculated in this work are systematically higher than those 
in A13. This is as expected, as the fainter DDP magnitude 
limit will include more galaxies in the calculation, resulting 
in higher densities. 

4.2 Sample Environments 

There are three key points to be investigated when compar¬ 
ing the environmental densities of ETG in the H-ATLAS 
sample with those from the HeViCS sample. Firstly, how do 


the respective sub-mm detected samples differ in environ¬ 
ment and is there any overlap? Next this study is extended 
to all ETG in A13 compared to all ETG in S13. Finally, 
it is also of interest whether the sub-mm detected versus 
non-detected ETG in these respective samples vary in en¬ 
vironment between themselves and if so, what the sense of 
this variation is. 

Trends of these distributions of densities are investi¬ 
gated in Figs. [4] and [5] A KS test of the sub-mm detected 
samples’ densities in Fig. ma) reveals a probability of only 
2x 10“^® of the two distributions being drawn from the same 
parent distribution. H-ATLAS and HeViGS ETG clearly re¬ 
side in very different environments, where HeViCS ETG are 
dominated by the dense cluster environment; the H-ATLAS 
ETG on the other hand mostly occupy sparse and non¬ 
cluster environments. 

Examination of the samples including those ETG with¬ 
out sub-mm detections in Fig. [3Kb) reveals the probability 
of ETG residing in the same environments is also zero, yet 
there is a modest overlap in the environments for the two 
samples between 20< Sgai <100 gals Mpc“^. This over¬ 
lap can mostly be associated with those ETG which are 
not detected at sub-mm wavelengths and is explored fur¬ 
ther in Fig. [S] Note that such an overlap is not apparent in 
Fig. Hi; a) for the sub-mm detected samples. The GAMA sur¬ 
vey as a whole is deep and wide enough to sample a broad 
range of galaxy environments, from isolate d held galaxies, to 
pairs , and both small and large groups ('e.g. lRobotham et al.l 
I 2 OIIII . However, it does not sample well the densest regions 
of the Universe as found in large clusters, since these are 
very rare environments. This can be seen in Fig.HKb), which 
shows that the GAMA galaxies in the three equatorial fields 
sampled in A13 do not extend up to the densities found in 
the Virgo cluster. Thus this ETG study is comparing and 
contrasting largely different environments. 

Fig.Oexplores the trend of surface density between sub- 
mm detected and undetected ETG for the respective sam¬ 
ples. KS tests for both sets of distributions indicate that 
there is a significant difference between the H-ATLAS dis¬ 
tributions (< 1% probability of them being the same), but 
that the HeViCS distributions are very similar. This indi¬ 
cates no environment density preference within the Virgo 
ETG (for the subsample of data used here). For H-ATLAS 
(Fig. HKa)), the sub-mm detected ETG have lower surface 
densities with respect to those of the undetected ETG - 
also shown in fig. 10 of A13. In Fig. HKb) for the HeViCS 
sample, the opposite effect is suggested. This latter result 
was noticed by S13 (and contrasted with the result for HI 
detections and non-detections), when fainter galaxies were 
included. The tendency found in HeViCS for dusty ETG 
to occupy the densest regions is consistent with the higher 
dust detection rate found for HRS ETG inside the Virgo 
Cluster than for t hose outside, particularly for lenticulars 
llSmith et al.ll2of3 l. 

This particular difference may be attributed to the fact 
that the two samples are environmentally very different (as 
explicitly shown in Fig. HI), with the H-ATLAS sample oc¬ 
cupying sparse environments and the HeViCS sample occu¬ 
pying a high density environment. However, given that both 
strangulation and ram pressure stripping in dense environ¬ 
ments are known to typically remove the ISM from galax- 
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Figure 4. Left panel: normalised distributions of sub-mm detected ETG in H-ATLAS (orange histogram) and HeViCS (black histogram) 
samples. Right panel: distributions of the combination of sub-mm detected and undetected ETG in H-ATLAS (orange) and HeViCS 
(black). 




Figure 5. Left panel: normalised distributions of sub-mm detected (orange histogram) and sub-mm undetected (blue histogram) ETG in 
H-ATLAS. Right panel: distributions of sub-mm detected (orange histogram) and sub-mm undetected (blue histogram) ETG in HeViCS. 
KS probabilities for whether the presented samples are drawn from the same distribution are shown in the top left of the plots. 


ie^, it is expected that sub-mm detected galaxies would be 
in lower density regions than undetected galaxies. The fact 
that the Virgo ETG are not exhibiting this behaviour indi¬ 
cates some other processes governing the presence of dust 
within these systems. Attention could be drawn to the case 
of M86: a Virgo elliptical which appears to have acquired its 
ISM via the strippi n g of gas and dust f rom a nearby spiral 
llGomez et al.ll20ig ). Smith et al.l ll2012l) also suggested that 
all their ETG acquired their dust through mergers. Gon- 
versely, b ased on the l a ck of e vidence for externally acquired 
material, iDavis et al.l ll2013l) argue against accretion as a 
general mechanism for gaining gas and dust in Virgo ETG. 

An additional effect which may be contributing to this 
difference is the ability of HeViGS to detect dust to lower lev¬ 
els than H-ATLAS: S13 amongst others found Virgo ETG 


® For example, dust stripping has been o bserved to be ongoing 
in the Virgo Gluster llGortese et al. I l2010all^ . 


with dust masses as low as 10® M©, and these lower dust 
masses appear to be quite common in the Virgo Cluster. 
Therefore by definition the dusty Virgo ETG are differ¬ 
ent to those being found by H-ATLAS. We observe that 
in the Virgo cluster there are no ETG with a 250/tm lu¬ 
minosity above 2.77x10^® W Hz”*^, which is equivalent to 
the threshold luminosity of H-ATLAS at the redshift upper 
limit (z=0.06). This is unexpected given that the samples 
are matched in optical luminosity (Mr <-17.4), however it 
does explain the differences in dust masses currently being 
observed. Therefore this difference in environments may be 
a cause of the differences in dust levels in these ETG (dust 
is destroyed in denser environments). 

Another possible cause for this difference that should 
be considered is the morphological classification of the ETG. 
HeViCS ETG have high enough optical resolution that they 
can be definitively categorised into their separate morpholo¬ 
gies. Given that H-ATLAS ETG lie at higher redshifts, their 
associated classifications cannot be assigned the same level 
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of acc uracy as the HeViCS ETG and in fact, iKelvin et alJ 
(|2014| group SO and Sa galaxies together in their clas¬ 
sifications. Since specihc dust ratio of galaxies systemati- 
call y increases when mov i ng from early- to l ate-type galax¬ 
ies dCortese et al.l l2012bl : ISmith et al. I 2 OI 2 I I , it is possible 
that a change in the threshold between ETG and late-type 
classification can skew the results. It is also well known 
that earlier-ty pe galaxies prefer denser environments (e.g. 
IPressletillOSOl ). Therefore any spurious LTG which may ex¬ 
ist in the H-ATLAS sample are likely to have both high 
dust-to-stellar mass ratio and sparser environments, thereby 
skewing the sample in the direction being seen. 

This last possibility can be investigated further by es¬ 
timating possible contamination le vels for the H-ATL AS 
sample. Original classifications from iKelvin et all ll20l4 l re¬ 
vealed 999 ETG, of which 285 were only agreed on by two 
of the three classifiers. Assuming that each classifier has ap¬ 
proximately equal weight, then these two-way agreements 
are estimated to be correct ~2/3rds of the time, leaving ~95 
incorrect classifications (LTGs misclassified as ETGs). Addi¬ 
tional criteria in A13 for creating the ETG samples included 
size and flattening (which removed 32 galaxies from the H- 
ATLAS sample, of which some may be LTG), plus evidence 
of spiral structure (which removed a further 22 galaxies from 
the H-ATLAS sample, which are all LTG). This reduces the 
potential number of contaminants to about 41 to 73 LTG. 
This is only 5 to 9 % of our sample of 771 ETG, whereas 
we detect 29% of our ETG sample in the submm. Therefore 
our observations are not consistent with resulting from LTG 
contamination alone. These values are estimates. This ques¬ 
tion of contamination of the H-ATLAS ETG sample will be 
better addressed in future through the use of deeper, sharper 
images fro m the VISTA and VST surveys covering GAMA 
areas (e.g. [Sutherland et al.l ((201^). 


5 MULTI-WAVELENGTH SED FITS 


PHYS is the physically realistic imposed energy balance be¬ 
tween absorbed starlight and dust emission. 

Emission from dust grains is modelled using a se¬ 
ries of modBB functions and a fixed template for poly- 
cy clic aromatic hydrocar bon (PAH) features, as described 
in Ida Cunha et ^ (l2008lL The free parameters in this dust 
model are the relative contributions of each component to 
the total IR emission and the temperatures of the warm 
circumstellar dust and cold interstellar dust We 

adopt an expanded version of MAGPHYS (da Cunha, priv. 
comm.) in the sense that the temperature ranges for warm 
and cold dnst are broadened to 30 K < < 70 K and 10 

K < < 30 K, respectively. This allows for a wider pos¬ 

sible range of temperatures found in some systems. It results 
in longer computation times, but permits better sampling of 
cold and low star-forming environments. 

A vast library is constructed by randomly drawing pa¬ 
rameter sets from the above model and constructing tem¬ 
plate SEDs with these sets. The expanded version of MAG¬ 
PHYS that we use here has 50k optical and 75k infrared 
templates in the library. For the dust an absorption coeffi¬ 
cient of K35 o= 4.54 cm^g“^, with /3 = 2 is assumed.The ob¬ 
servational SED is then modelled b y comparing a li b rary o f 
stochastic models (as described in Ida Cunha et al.l (l2008l H 
to the observed data and weighing the output parameters 
with the corresponding yf, constructing probability distri¬ 
bution functions (PDFs). Depending on the data coverage of 
the SED, some parameters are more accurately constrained 
than others. In this paper, we limit ourselves to the param¬ 
eters listed below: 

Cold and warm dust are responsible for the large part 
of the total dust mass in galaxies (M^^*^ and M(^*^ re¬ 
spectively). Warm dust in birth clouds (M^'^) can also con¬ 
tribute. The contributions of hot dust and PAHs are ac¬ 
counted for in a multiplicative factor of 1.1. Thus the total 
dust mass is the sum of these components: 


The photometric data of both the HeViCS and H-ATLAS 
sample cover a wide wavelength range. Intrinsic informa¬ 
tion is encoded in the spectral energy distribution (SED) 
of each galaxy. Certain wavelengths are directly linked to a 
single component (e.g. NIR emission traces the old stellar 
population), while others are ambiguous (e.g. optical light 
is influenced by both stars and dust). It is therefore useful 
to treat the multi-wavelength information of a galaxy at the 
same time using a complete model. 

5.1 MAGPHYS 

MAGPHYS - Multi-wavelen gth A nalysis of Galaxy Physical 
Properties (Ida Cunha et al.l 2008ll is a Bayesian fitting code 
which is able to model the UV to snbmm SED of galaxies. 
The program relies on a multi-component galaxy model to 
predict the flux in eac h wavelen^h. S tarting from an ini¬ 
tial mass function from IChabrieil I 2 OO 3 I I . stellar components 
are evolved in time using the stellar pop ulation synthesis 
(SPS) model of iBruzual fc Char loti (l2()03l i. The interaction 
of starlight with diffuse interstellar and star -forming region 
dust i s calculated from the two-component ICharlot fc Falll 
(l200(tl extinction model. One of the key points of MAG- 


Cold dust in the diffuse ISM has an equilibrium temperature 
represented by 

Accordingly, the equilibrium temperature of warm dust in 
birth clouds is represented by . 

The total amount of infrared light emitted by dust grains is 
parametrised in the total dust luminosity Ldust- 
The total stellar mass M* as derived from the SPS models. 
The star formation rate (SFR) is an average of the mass of 
stars formed per year, during the past 100 Myr. The under¬ 
lying star formation law is an exponentially declining SFR 
starting from the birth of the galaxy. Throughout the life¬ 
time of the galaxy, there is a random chance of starbursts 
taking place. 

The ratio of the SFR and M« is then called the specific star 
formation rate (sSFR). 

The time at which the galaxy is formed, Tform, is defined as 
the age of the oldest stellar population. 

The time at which the last starburst ended, Tiastb. 

5.2 Data coverage of the SED 

Good coverage of the panchromatic SED is desirable to en¬ 
sure reliable modelling. The GAMA/H-ATLAS dataset com- 
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Figure 6. MAGPHYS rest-frame SED fit to H-ATLAS elliptical galaxy 298980. The SED is fit to observed photometry in the FUV, 
NUV, SDSS ugriz, UKIDSS YJHK, WISE Wl, W2, W3 and W4, PACS 100, 160pm and SPIRE 250, 350 and 500pm wavebands (black 
points). The red line shows the overall attenuated model fit, whilst the blue line shows the unattenuated optical model. Below are the 
residuals from the fitted points, and below these are the likelihood probability functions for parameters of this elliptical galaxy. 


prises of self-consistent photometry based on a standard 
format. Therefore the H-ATLAS ETG sample has the fol¬ 
lowing data coverage of the SED: GALEX FUV and NUV, 
SDSS ugriz, UKIDSS YJHK, WISE W1-W4, PACS 100 and 
160pm and SPIRE 250, 350 and 500pm. Note that some 
galaxies are missing data as detailed below: 

GALEX FUV: 38 galaxies (17%) 

GALEX NUV: 36 galaxies (16%) 

WISE W1-W3: 5 galaxies (3%) 

WISE W4: 90 galaxies (41%) 

PACS 100pm: 19 galaxies (9%) 

PACS 160pm: 26 galaxies (12%) 

Missing galaxies in these wavebands are due to an inability 
to match the optical source to a counterpart in the specific 
waveband (i.e. no detection). The 22pm WISE W4 band 
in particular suffers from a low detection rate due to the 
waveband’s low signal-to-noise. Data in other wavebands is 
complete for all the galaxies. 

In the case of the HeViCS ETG, the UV part of the 
spectrum is covered by fluxes from the GALEX catalogue 
except for t he 12 galaxies which a re also in HRS. For those, 
fluxes from ICortese et al.l (l2012bl) were used. Optical and 
NIR filt ers include UBV fc JHK from the GOLDMINE 
archive (iGavazzi et al.llioosl l. 


WISE dWright et al.ll^ld) detected almost all galax¬ 
ies in the HeViCS sample during its all-sky survey. Unfor¬ 
tunately archival fluxes caused inconsistencies in the SED 
fits. The MIR is an ambiguous regime, with both emission 
from old stars and hot dust. Sufficient sampling and correct 
flux determination are vital to disentangle both components 
in an SED fit. We therefore chose to redo the flux mea¬ 
surements on the archival WISE images and minimise the 
contamination of foreground stars. Appendix 1X1 gives more 
detail on the measurements. 

Herschel observations at 100 — 500 pm were taken from 
S13 and complete our dataset. 

5.3 SED Results 

MAGPHYS was used to fit energy balance models to each of 
the HeViGS and H-ATLAS ETG, as described above. Figs. 
[6] and [7] show example MAGPHYS fits to one of the H- 
ATLAS and HeViGS ETG respectively, with the resultant 
PDFs for a variety of fit parameters shown in the lower 
panels. MAGPHYS cannot include a synchrotron compo¬ 
nent and is therefore unable to correctly fit the four dusty 
HeViGS ETG with such a component, and therefore these 
are excluded from the multi-wavelength analysis. 

In order to gain some insight on the goodness-of-fit 
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Figure 7. MAGPHYS rest-frame SED fit to HeViCS elliptical galaxy VCC 408. The SED is fit to observed photometry in the FUV, 
NUV, UBV, 2MASS YJHK, WISE Wl, W2, W3 and W4, PACS 100, 160/rm and SPIRE 250, 350 and 500/rm wavebands (black points). 
The red line shows the overall attenuated model fit, whilst the blue line shows the unattenuated optical model. Below are the residuals 
from the fitted points, and below these are the likelihood probability functions for parameters of this elliptical galaxy. 


for each galaxy, the 29 (33 minus 4 synchrotron galaxies) 
HeViCS fits were visually inspected and assigned a flag for 
‘good’ or ‘poor’ fit. Four galaxies were assigned ‘poor fit’ sta¬ 
tus - each of these fits had an associated reduced valueQ 
greater then 4. This was then chosen as the criterion to as¬ 
sess whether the H-ATLAS fits were ‘good’ or ‘poor’. Eleven 
H-ATLAS systems were found to have ‘poor fits’. These fif¬ 
teen galaxies from the two samples are also highlighted in 
future plots to separate them from the rest of the sample. 
Additionally, all ETG with ‘poor fits’ or a synchrotron radio 
component are excluded in any further statistical analysis in 
this section. 


5.3.1 Contrasting Derived Parameters 

As described in ll3.1l and 113.21 A13 and S13 fit single modBBs 
to their FIR/sub-mm data to obtain dust masses and tem¬ 
peratures for their ETG. These dust masses are used to 
evaluate the dust-to-stellar mass ratio which is plotted as 
a function of stellar mass in the left panel of Fig. [S] The 


® The MAGPHYS 1® a constraint on the best fitting theoret¬ 
ical template SED and hence the most likely fit. 


stellar masses are from AI3 and S13. The effects of detec¬ 
tion limits are illustrated. The red dashed lines in Fig. [8] 
represent the range of dust mass detection limits at the 250 
micron flux limit of 33.5 mjy for the H-ATLAS sample, with 
the lower line for nearby (z=0.013), high temperature (30K) 
detection limits (as shown in A13) and the upper line for 
far (z=0.06), low temperature (I5K) detection limits. This 
includes the typical temperature range of most H-ATLAS 
detections. Below the lower line no detections are expected 
and above the upper line all sources in this temperature 
and redshift range are detectable. Between the two lines de¬ 
tectability depends on the distance and dust temperature 
of sources, with decreasing likelihood of detection going to¬ 
wards the lower line. The blue dot-dashed lines represent the 
corresponding limits for the HeViGS 250 micron flux limit 
of 25.4 mJy, the distance range of HeViCS sources (17 to 32 
Mpc), and considering the effects of different dust temper¬ 
atures for the various morphological types from section 4.1 
of S13. The specific dust mass detection limits are lower for 
the HeViGS sample, which is nearby, mostly at a distance of 
17 Mpc. Fig.[8]shows a key difference in the normalised dust 
levels of the two ETG samples. It indicates that the HeViCS 
ETG have less dust, by a factor of 10 or more, than the H- 
ATLAS ETG. This difference is due to the much lower dust 
detection limit of HeViCS, whose galaxies are approximately 
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ten times closer than the H-ATLAS galaxies. However, it re¬ 
mains to be understood why there are no HeViCS ETG as 
dusty as the dusty H-ATLAS ETG, at fixed stellar mass. 

A two sample KS test can be applied to the left 
plot in Fig. [8] by considering a diagonal line parallel to 
the limits and above which there are no HeViCS detec¬ 
tions. This line goes through VCG1535 (at logi0(M*)=10.94 
and logio(Md/M*)=-4.28 in Fig. [8] left plot). 162 H- 
ATLAS/GAMA ETG lie above this line, out of 771 GAMA 
ETG with Mr < — 17.4mag. For the Virgo sample there 
are 123 (including all detected and undetected) ETG with 
Mr < — 17.4mag in the HeViCS area, which all lie some¬ 
where below the diagonal line. Therefore the KS statistic is 
0.21, which for these optical sample sizes has a probability 
< 0.1% of them being the same. Thus the two samples differ 
significantly from each other in Fig.|5] More studies of ETG 
in other clusters would help to verify if this is a difference 
due to environmental density. 

Normalised dust mass calculated from MAGPHYS is 
shown plotted against stellar mass in the right panel of 
Fig. [8l where the stellar masses are also from MAGPHYS. 
Both ModBB and MAGPHYS plots are shown side-by-side 
for ease of comparison. Qualitatively there is little difference 
between these distributions. There is a negative trend found 
for this relationship for both these samples, partly driven by 
selection effects in the lower left corner. 

Dust properties from MAGPHYS fits are examined and 
contrasted to the modBB solutions in Figs. [Oja) and [9Kb). 
The MAGPHYS most likeljQ cold dust temperature and 
overall dust mass is shown on the x-axis, and modBB so¬ 
lutions are shown on the y-axis for both samples of ETG. 
As previously stated, poor fits are shown as different sym¬ 
bols to separate them from the good fits. 

Fig.EKa) shows that dust temperatures are similar, ex¬ 
cept for a few outliers that vary more significantly between 
the fitting methods. These tend to be for galaxies without 
PACS detections. The scatter observed here is in the sense 
that MAGPHYS assigns higher temperatures to the likely 
cold dust grain distributions. This difference may result from 
the fact that MAGPHYS fits multiple temperature compo¬ 
nents, unlike the single component ModBB. More MIR data 
coverage would help to constrain the dust temperatures in 
these outliers. However these differences in temperatures do 
not necessarily cause similar scatter in other derived param¬ 
eters. 

MAGPHYS dust mass shows a better correspondence 
with dust masses derived from modBB fitting (see Fig.[9Kb)). 
There is a slight offset for some galaxies in both samples from 
the x=y plane - typically the modBB fitting appears to give 
higher dust masses. This is because MAGPHYS assumes a 
higher on average dust temperature than modBB fitting, 
thereby resulting in lower dust masses. 

These tests reveal that the dust mass parameter is well 
described by both a modBB and MAGPHYS, as the mean 
results for the two samples do not change substantially, nor 
do the KS-test probabilities differ between the two meth¬ 
ods. However, a small difference is found between the stel¬ 
lar masses derived by GAMA/HeViGS teams and those de- 


^ The most likely value of a parameter is chosen as the median 
of the probability distribution functions output by MAGPHYS. 



Figure 9. Comparison of derived cold dust temperature (top), 
dust mass (center) and stellar mass (bottom), where the y-axis 
represents values from modBB fitting and the x-axis represents 
MAGPHYS fitting results. For the stellar mass comparison, the 
y-axis represents stellar mass values from A13 and S13. The H- 
ATLAS sample is plotted in red and HeViCS ETG are plotted in 
blue. Galaxies with poor fits (x^ >4.0) are shown as open squares. 
The solid line is a one-to-one line and error bars represent the 
mean la errors in either direction for both samples. Inset plots 
show the difference offset from zero for each parameter. 
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Figure 8. Dust-to-stellar mass ratio plotted as a function of stellar mass calculated using ModBBs (left panel) and MAGPHYS (right 
panel). HeViCS (blue points) and H-ATLAS (red points) samples are shown in both plots, and galaxies are subdivided into E (red squares 
for H-ATLAS and blue squares for HeViCS) and SO (red dots for H-ATLAS and blue dots for HeViCS) classifications. In the MAGPHYS 
plot, galaxies with poor fits (PF) are encircled crosses in the samples’ respective colours, while all other ponts are good fits (GF). Error 
bars in the left panel give the mean overall uncertainty on the points from ModBB fits, in the same colours as their respective samples. 
In the left plot diagonal lines represent the range of dust mass detection limits for the samples in their respective colours (see Section 
5.3.1 for details). These same lines are included as a guide only in the right plot. Error bars in the right panel give Itr to each side of the 
PDF. 


rived by MAGPHYS. Fig. 9(c) shows that stellar masses ob¬ 
tained with MAGPHYS are on average slightly smaller (by 
~ 0.1 dex) th an those obt a ined previously, particularly at 
small masses. iTavlor et al.l 1 I 2 OIII ) showed that inclnsion of 
UKIDSS fluxes slightly biases stellar mass measurements in 
this way. This is a small effect for the cnrrent applications. 
No systematic difference between HeViGS and H-ATLAS 
ETG is apparent in Fig. 9(c), which demonstrates that the 
previous stellar masses were also obtained consistently for 
the two samples. 


We compare the relationship between derived dust and 
stellar mass for the two samples in Fig. IIUI Trends are 
found for both H-ATLAS and HeViGS ETG, with corre¬ 
lation coefficients of rp=0.42 and 0.58 respectively, both 
with < 1% probability of chance occurrence. The differ¬ 
ence in dust mass between the samples is highlighted by 
the ranges exhibited: 4.1<log(M(j/M0)<6.7 for HeViGS and 
5.5<log(Md/MQ)<7.9 for H-ATLAS. Results from stacked 
data for red galaxies in (iBourne et al.lFoi2l . their fig. 16) 
show a similar trend to that of the H-ATLAS ETG shown in 
Fig. IIUI Following a similar KS test as previously, but using 
a line through the upper edge of the HeViCS detections in 
Fig. HOI we And 197 H-ATLAS detections above the line, out 
of 771 GAMA ETG. This gives a KS statistic of 0.25, which 
for these optical sample sizes gives a probability < 0.001% 
of the Virgo ETG being drawn from the same parent pop¬ 
ulation as the GAMA ETG. In other words the proportion 
of dusty ETG is much higher in the GAMA sample than in 
the Virgo sample. 



Figure 10. MAGHPYS dust mass plotted as a function of stellar 
mass. HeViCS (blue points) and H-ATLAS (red points) samples 
are shown, where galaxies are subdivided into E (red squares for 
H-ATLAS and blue squares for HeViCS) and SO (red dots for H- 
ATLAS and blue dots for HeViCS) classifications. Galaxies with 
poor fits are encircled crosses in the samples’ respective colours. 
Error bars give Icr to each side of the PDF, in the same colours 
as their respective samples. 


6 DISCUSSION 

6.1 Mock Catalogues and Detection Limits 

We further test how the two samples of ETG compare in 
the left plot of Fig. [8] by carrying out Monte Carlo simula¬ 
tions for the Virgo galaxies, to account for the sample selec- 
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tion effects of the GAMA galaxies. 500 mock catalogues, 
each with 771 galaxies (as in the GAMA ETG sample), 
are generated by randomly sampling the Virgo ETG (for 
their dust mass, stellar mass and submm fluxes) and placing 
them at distances randomly selected from the GAMA ETG 
sample. The observed Virgo ETG submm fluxes (or non¬ 
detections) in 250 and 350 micron Herschel bands are then 
transformed to fluxes at the new distances, taking into ac¬ 
count K-corrections (as in equation 2 of Dunne et al. 2011). 
Thus fluxes significantly decrease due to the larger distances 
and slightly increase due to K-corrections in the submm. 

For each mock catalogue the number of expected 
submm detections above the H-ATLAS limits (5 sigma at 
250/rm and 3 sigma at 350/im, as in A13) is calculated and 
expressed as a percentage of the 771 galaxies. The predicted 
mean and standard deviation is then 0.92 ± 0.36% detec¬ 
tions, with standard deviation calculated from the spread 
of results amongst the 500 mock catalogues. For GAMA 
ETG 188 out of 771 galaxies were actually detected (above 
those H-ATLAS limits, as in Fig. [8] left). This corresponds 
to 24.4 ± 2.0%, taking into account Poisson noise in these 
numbers. 

Therefore the difference between H-ATLAS ETG de¬ 
tections and the mock HeViGS detections at H-ATLAS dis¬ 
tances, is greater than 11 sigma. This corresponds to a neg¬ 
ligible probability (< 10“^®) of the HeViGS sample being 
consistent with the H-ATLAS sample, in terms of their in¬ 
trinsic submm properties. Thus the differences in specific 
dust distributions seen in Fig. [8] is not explained by differ¬ 
ent dust mass detection limits in the two samples. 

6.2 Environments 

As it has been shown that MAGPHYS successfully repro¬ 
duces the modBB results for dust mass, we adopt hereafter 
the parameters (Ti, Ma and M*) derived with MAGPHYS. 
This has the advantage of ensuring uniformity of approach 
and also of providing SF parameters. However, this means 
that we must exclude from the subsequent analysis the 4 
HeViGS ETG with synchrotron radiation and the MAG¬ 
PHYS poor fits (4 in HeViGS and 11 in H-ATLAS). We note 
that the multi-waveband photometric datasets used in the 
MAGPHYS fits are different for the HeViGS and H-ATLAS 
samples. We begin by examining the relationship between 
environment and dust-to-stellar mass ratio in Fig. 1111 In 
spite of the large scatter in this plot, there is a weak overall 
log-log anticorrelation (rp=-0.3), which is influenced by the 
different sensitivity limits of the two samples. For H-ATLAS 
data alone the correlation coefficient decreases to rp=-0.2, 
which is not a significant result, therefore the anticorrelation 
is not clear within each individual sample. 

For the H-ATLAS sample, which is less sensitive, we 
have attempted a stacking analy sis to probe th e unde tected 
ETG (following the methods of I Bourne et al.l (l2012l lL ac¬ 
counting for blending of nearby sources and assuming the 
flux is spread o ut over the optical extent of each source 
dHill et al.ll201ll L By stacking at the position of the unde¬ 
tected ETG in the 100 - 500 /rm maps we derive the median 
value of the flux density in the stack in order to avoid bias 
from outliers, and we estimate the la error on the median 
from the distribution of values in the stack. From these me¬ 
dian stacked flux densities we find a median dust mass of 



Figure 11. Dust-to-stellar mass ratio derived from MAGPHYS 
fits plotted against environment surface density for the two ETG 
samples. H-ATLAS ETG are plotted in red (Es are open circles 
and SOs are filled circles) and HeViGS ETG are plotted in blue 
(Es are open squares and SOs are filled squares). Poor SED fits are 
marked in encircled crosses (in red or blue based on the sample). 
Arrows indicate the dust-to-stellar mass ratio for two H-ATLAS 
SOs at Sgai =0.001 gals Mpc“^. 


4.44x10®Mq, at a median surface density of Egai = 2.3 
gals Mpc“^ and median stellar mass of 2.0 x 10^°M©. 
This is at logio(M*) = 10.3 and logio(Mjj/M*)=-4.65, which 
lies amongst the dust detected SO galaxies in the Virgo clus¬ 
ter in Fig. |5] and well above the HeViGS detection limits. In 
Fig.HHthis specific dust mass is comparable to the Virgo de¬ 
tected ETG. However the 90 undetected ETG in the Virgo 
cluster will lie well below this point in specific dust mass, 
since they are below the dashed blue lines in Fig. [S] and the 
two samples cover similar stellar mass ranges. Gonversely, 
if the detected Virgo ETG were placed at the average dis¬ 
tance of the GAMA ETG sample, then only one Virgo ETG 
would be detected in Fig. [TT] (VGC1535). It is difficult to 
quantify the lower limit of the distribution in Fig. [TT] with¬ 
out further constraints, but it is clear that sensitivity limits 
are not causing the downward slope seen at the top of the 
distribution in Fig. 1111 

The substantial range of normalised dust mass dis¬ 
played by the Virgo ETG is a feature which needs to be 
subjected to more scrutiny, e.g. studying the dust detection 
rate and normalised dust mass in fast- and slow-rotators, fol¬ 
lowing a suggestion by S13. ATLAS®^ is an ongoing survey 
investigating the kinematic properties of a volume-limited 
sample of ETG including the Virgo cluster, finding that el¬ 
liptical galaxies tend to be either slow rotators or rotate 
faster than lenticulars (based on apparent specific angu¬ 
lar momentum). They find that non-rotati ng ETG tend to 
be fo und in highly overdense environments (iKrainovic et al.l 
[2^ - their results also indicate that in dense groups and 
clusters gas accretion is suppressed. Although twenty of the 
ETG in the H eViGS sample corresp ond with those studied 
by ATLAS®^ (lEmsellem et al.|l201l] ). it is inconclusive from 
such a comparison whether the rotation speed of each galaxy 
is related to its respective dust-to-stellar mass ratio, for the 
current sample. However it has been shown that a galaxy’s 
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stellar angular momentum and stella r mass are negatively 
correlated fe.g. lEmsellem et al.ll20l'ill . and normalised dust 
mass is negatively correlated with stellar mass (A13). There¬ 
fore it can be postulated that dust-to-stellar mass ratio is 
positively correlated with a galaxy’s rotation speed. How¬ 
ever, the reverse trend was found in Virgo (S13), probably 
because in a denser environment the giant slow rotators may 
be collectors of stripped dust. This suggestion might be in¬ 
vestigated in the future with high-resolution, kinematic ob¬ 
servations of the cold dust or gas in GAMA/H-ATLAS ETG. 

Given the wide range of environments inhabited by the 
H-ATLAS ETGs, and the large sample size, we can now 
split H-ATLAS ETG into two subsamples: H-ATLAS-LOW 
(Egai <1.25 gals Mpc-=) and H-ATLAS-HIGH (Sg^i >1.25 
gals Mpc“^jj. Note that such subsampling is not possible for 
HeViCS ETGs due to the small sample size. We examine 
average MAGPHYS dust-to-stellar mass ratio for the two 
subsamples, and subsequently compare with HeViCS dust- 
to-stellar mass. 

H-ATLAS-HIGH has a mean dust-to-stellar mass ra¬ 
tio of logio(Md/M*)=-3.35, indicating lower dust frac¬ 
tions on average in comparison to H-ATLAS-LOW (mean 
logio(Md/Mt)=-3.19). This may be influenced by the well 
known property that massive ETG tend to reside in denser 
environments. In fact H-ATLAS-HIGH has intermediate 
dust-to-stellar mass ratio between H-ATLAS-LOW and 
HeViCS (mean logio(Md/M*)=-4.55). However, a KS test 
indicates no significant difference in the distribution of this 
parameter between the two H-ATLAS subsamples. 

These results for dust parameters can be interpreted 
in different ways. The subtle decrease in relative dust mass 
with increasing environment may be a real effect, which is 
then extended to the dense region of the Virgo cluster, or 
it may be completely spurious. Our tests using mock cat¬ 
alogues above confirms the real difference between the two 
samples. More samples at intermediate environment densi¬ 
ties and in other cluster environments are needed to test the 
reality of this trend with environment. 

6.3 Star Formation Properties 

In A13 UV-optical colour was used as a proxy for SFR in a 
galaxy. However blue UV-optical colours can also be induced 
by the presence of a very old stellar population in the galaxy 
dGreggio fc Renzhiil Il990l : iHorch et alj 1 19921 : iBressan et al.l 
Il994l ~) and therefore further investigation of the SFRs in 
these galaxies is a necessity to confirm the results found thus 
far. This must be handled with care, as SFRs derived from 
MAGPHYS are also related to the UV emission; however 
the inclusion of longer wavelength information and energy 
balance in SFR calculations gives a better estimation than 
a simple proxy. Additionally, it should be noted that SFR 
estimates using FIR data are subject to large errors in galax¬ 
ies where dust heating is dominated by the diffuse radi ation 
field from an old stellar population llBendo et al.ll20l^ . 

The interpretation of SED fits to ETG data must ac¬ 
count for the potential contribution to UV light from old, 
evolved stars on the horizontal branch ('HBi. fBurstein et al.l 

® The environment cutoff has been chosen to create subsamples 
of equivalent size 


lll988ll showed that populations of stars older than ~10 Gyrs 
can cause an upturn in the UV flux, due to UV emission 
from HB stars. Younger populations do not have this com¬ 
ponent contributing to the integrated UV light. Populations 
younger than -^3 Gyrs again have excess UV emission but 
from the young, massive, main sequence stars. Thus all in¬ 
termediate age populations (-^3 to 10 Gyr) do not have sig¬ 
nificant UV upturns in their spectra and hence will show 
no sign of a UV excess that could be erroneously attributed 
to the presence of young stars. The GALEX NUV flux is 
less affected by the UV flux from old, evolved stars than the 
GALEX FUV band fsee lKavirai et aT]|2007l . their fig. 1). 

A concern with using MAGPHYS for the interpretation 
of parameters such as SFR, is that they are based on opti¬ 
cal, UV and FIR emission and are calibrated only for galax¬ 
ies where dust heating is primarily contributed to from the 
young and old stellar population. However, if the photons 
heating the dust come primarily from UV emission from an 
old stellar population, then the results obtained with MAG¬ 
PHYS are less reliable, because of the uncertainty in UV 
contribution from old stars. This is a known effect, which 
is part ly incorporated into the models of iBruzual fc Charlotl 
(l2003ll used in MAGPHYS, throu gh the (albei t uncertain) 
inclusion of some hot evolved stars. I^lim et aP (l2007fl (and 
references therein ) tested nearby elliptical g alaxies, also us¬ 
ing the models of iBruzual fc Gharlod (l2003l l and found that 
those models could account for UV light from old stellar 
populations. In addition, eyeballing the fits done here indi¬ 
cates that the UV upturn of the ETG is not a strong effect. 

We can further strengthen our argument that UV emis¬ 
sion from the old stellar population is not the driving mech¬ 
anism for the dust heating by examining the NUV-r colour, 
in relation to that expected from the UV upturn. UV con¬ 
tamination from old stars leads to colours that are still 
redder than NUV-r>5.0 (from the sample of iKavirai et al.l 
(I 2 OO 7 II : their fig. 11). Their ETG classified as old (ages>l 
Gyr from stellar population model fits) all have NUV-r>5.0. 
For example, NUV-r = 5.4 is t he colour of the strong UV- 
upturn galaxy NGG 4552 fsee lYi et al.l lioosl i. In the H- 
ATLAS/GAMA sample 141 out of 184 NUV detected ETG 
are bluer than NUV-r=5.0. This large fraction of blue NUV- 
r ETG in the A13 sample indicates that the blue colours 
cannot be explained by a UV upturn from old, evolved stars 
alone. Therefore the blueness and large scatter of their NUV- 
r colours (ranging from ~1.5 to 6.5) indicates that these 
colour and UV fluxes are dominated by different amounts 
of recent star formation in these ETG rather than by UV 
emission from old, evolved stars. However, we note that for 
up to 36% of the 220 dust detected ETG in A13 (includ¬ 
ing 43 with NUV-r>5.0 plus 36 NUV non-detections) the 
UV radiation and dust heating could be dominated by old 
stars, making their MAGPHYS SF rates overestimates. The 
fractional mass involved in the star formation does not need 
to be very large in order to strongly influen ce the NUV-r 
colours of stellar populations (see review by iKavira ][200i, 
their fig. 1). 

The 33 ETG in the HeViCS sample do not exhibit 
such blue colours , with the average NUV-r'~5.3. Also 
Idi SeregQ Alighi^ ll2013h has shown that dust-detected 
HeViCS ETG are not bluer in B-H colour than the unde¬ 
tected ones. These facts should be kept in mind when con¬ 
sidering the parameters extracted based on MAGPHYS fits 
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Figure 12. MAGPHYS derived specific star formation rate plot¬ 
ted against stellar mass for the two ETG samples. See Fig. 1111 
for symbols and labels. Green stars overplot points with NUV- 
r>5. Red and blue dashed lines show best linear fits to H-ATLAS 
and HeViCS respectively. Histograms of these distributions are 
plotted on the right. 


to these galaxies, and may in fact result in star formation 
rate overestimates for Virgo Cluster ETG. 

We examine the sSFR (defined in Section [57I|l as derived 
by MAGPHYS below. Fig. [12] shows sSFR plotted against 
stellar mass, where a similar trend to Fig. [ 8 ] emerges in the 
form of an anticorrelation between sSFR and stellar mass of 
the ETG for the two samples. Regression lines are fit to the 
two samples, revealing correlation coefficients of xp =-0.572 
and -0.733 for H-ATLAS and HeViGS samples respectively. 
These are both significant with much less than 1% prob¬ 
ability of occurring by chance. As expected, since derived 
dust mass and SFR are both influenced by the FIR flux, 
this corresponds to what was observed for dust-to-stellar 
mass: H-ATLAS ETG show a weaker correlation for dust- 
to-stellar mass and sSFR against stellar mass in comparison 
to HeViCS. Those galaxies with NUV-r>5.0 are indicated 
by green stars in Fig. 1121 This gives an indication of which 
points are most likely to be upper limits rather than full SFR 
measurements in this plot and most of these occur at higher 
stellar masses. Future work is required to more accurately 
constrain the low specific star formation rates occurring in 
massive ETG. 

The side panel of Fig.[T2]shows the distribution of galax¬ 
ies in sSFR space for the two samples. This quite clearly 
identifies the higher-on-average sSFR for H-ATLAS ETG 
compared with HeViGS ETG, further strengthening the 
point that H-ATLAS ETG are not only dusty, but actively 
star forming. The most star forming end of this distribution 
displays sSFRs which exceed that of our Milky Way Galaxy. 

This anticorrelation between sSFR and stellar mass has 
previously been o bserved in both t he local and medium- 
redshift Unive rse l|Salim et al.llioo^ Somerville et al.ll200^ : 
iFirmani et al ] l20inl L The study bv I Salim et al.l |2003) took 
observations of 50,000 SDSS galaxies with a range of mor¬ 
phologies and stellar masses and, after measuring their sS¬ 
FRs using synthetic population models including dust at¬ 
tenuation, constrained this relation to purely star-forming 


galaxies as: 

log(sS'Fi?) =—0.35 log M, — 6.33. (3) 

Those star-forming galaxies covered log(M«/M 0 ) ~8.4 to 
11.3 and log{sSFR) ~- 8.6 to -10.9. In contrast with this, 
and over a similar range in stellar mass, the HeViCS ETG 
display a slightly steeper slope with: 

log(sS'FR) = -0.59log M, - 6.39 (NUV-r = all) (4) 

log{sSFR) = -0.39 log M* - 8.02 (NUV-r < 5.0) (5) 

and the H-ATLAS ETG produce the steepest gradient of all 
with: 

log(sSFR) = -1.37log M* -b 3.40 (NUV-r = all) ( 6 ) 

log(s5'Fi?) =—0.90 log M* — 1.20 (NUV—r < 5.0) (7) 

These differences may be attributed to both the larger sam¬ 
ple size and the different galaxy types (presumably dom- 
inated by later-ty pe galaxies) making up the relation in 
ISalim et al.l (120071 ). For the HeViCS and H-ATLAS samples 
we give above the relations both with and without the green 
points shown in Fig. 1121 fNUV-r>5.0i. since those points in¬ 
dicate the least certain associations with ongoin SF. Irre¬ 
spective of whether those points are included or not, the 
ETG used in our work span a large range of sSFRs and 
the relation’s gradient is likely to be steepened by the ex¬ 
treme sSFR values. Such a steep relation for the H-ATLAS 
ETG, coupled with the extreme levels of dust content for 
lower stellar mass H-ATLAS ETG, is consistent with sSFR 
downsizing where lower mass ETG harbour star formation 
in even the local Universe. These low mass ETG also oc¬ 
cupy the sparsest environments (<1 galaxy Mpc“^), further 
strengthening this downsizing theory and in accordance with 
previous results shown for galaxies in the local Universe 
llCassata et al.l l2007l : ICooner et al.l 120071 : IWiiesinghe et al.l 

I 2 OI 2 II . 

What we are seeing in the H-ATLAS E TG sample are 
galaxi es like the r ejuvenated on e s prop osed in iThomas et al.l 
l201fll l. (see also iVoung et al.l (120141 ')!. These rejuvenated 
ETG have some recent star formation as well as old stars 
and are increa s ingly common with decreasing galaxy mass. 
iThomas et al.l (l201f)l ! postulate, from trends seen in SDSS 
data, with environment and chemistry, that these cases, with 
contributions from relatively recent star formation, are deci¬ 
sively influenced by environment. They support their claim 
of rejuvination with observations that show that the more 
recent star formation is less enhanced in [a/Fe], and has 
thus had time to build up more iron from type-la super¬ 
nova, for the composition of the later starburst, in contrast 
to the older stars. This trend with environment is similar 
to the results that we find here. Their fig. 8 shows that this 
is particularly significant for low mass ETGs, as we find in 

Fig- El 

Models for downsizing predict that lower mass galax- 
ies have more exte nded star formation histories (e.g. 
IPe Lucia et al.l l2006h . If the H-ATLAS sample is divided 
into low mass (M* < 1 O^°'®M 0 ) and high mass systems 
(M* >1O'^° ®M 0 ), it is interesting to note that the lower 
mass systems exhibit higher average sS FRs over the full red - 
shift range (consistent with results bv IFirmani et alll 201 f)l !. 
Additionally the higher mass systems show the most similar 
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sSFR to the average sSFR of the HeViCS ETG, indicat¬ 
ing greater similarity between these systems than the lower 
mass ETG have with either of these groupings. 

6.4 Age Properties 

Ideally, an exploration of the ages of these ETG would begin 
by using spatially resolved, population synthesis modelling 
for these systems as a whole. In the case of the H-ATLAS 
ETG this has not been done yet because current imaging 
does not have good enough resolution. Therefore we choose 
to run a pilot study on the ages of these galaxies using MAG- 
PHYS results. Caution must be applied to the use of these 
results, as they are fully dependent on SPS code used to com¬ 
pute the short-wavelength light produced by stars, which is 
likewise dependent on the model’s choice of metallicity, ini¬ 
tial mass function (IMF) an d star formation his t ory. I n this 
case, the SPS code is that of iBruzual fc Chariot! l|2003l 'l. and 
they adopt a range of expo nentially declin ing SFFls plus 
bursts and a Chabrier IMF (|Ghabrierll2nn,^ . Here we will 
only consider relative numbers, and not absolute ages. 

We choose to examine two forms of galactic age: the for¬ 
mation timescale (Tform), which is defined as the age of the 
oldest stars in the galaxy and is a good representative of the 
age of the galaxy, and the time the last burst of star forma¬ 
tion ended (Tiastb). The distributions of both these parame¬ 
ters and means thereof are compared in Fig. 1131 Probability 
results from KS testing the distributions are also included in 
these figures. The left panel indicates that the two samples 
have significantly different formation timescales. However, 
the results for Tiastb indicate that there is no difference (<1 
per cent probability of a difference) between the two sets of 
galaxies. This is an interesting result as it is the first point at 
which any similarity between the parameters of the two sets 
of ETG has been found, although Tastb is not always well de¬ 
termined. The green histograms plotted for Tastb in the right 
plot in Fig.[T3]illustrate where the least certain meaurements 
are, since these have NUV-r>5.0 colours. These timescales 
are provisional estimates, compared here in a relative way. 
Further detailed study on the galactic stellar populations 
is required to determine whether these results are real or 
simply a result of the assumptions made in the SPS fitting. 


6.5 Additional Considerations 

Based on the results discussed above, a clear conclusion 
about the two samples is reaffirmed: that the ETG in each 
of the samples have differing dust properties, with HeViCS 
ETG demonstrating consistently low dust levels, whilst the 
H-ATLAS ETG have significantly higher dust levels which 
bridge the gap between HeViGS ETG and late-type spirals. 
The dust-to-stellar mass ratio is shown to be strongly driven 
by the stellar mass of the galaxy, particularly for the HeViCS 
systems. This work has made uniform the calculations of 
stellar and dust mass for two samples; this has served to 
strengthen the result that there is no overlap between dust- 
to-stellar mass ratio for fixed stellar mass for ETG from 
these two samples, covering different environments. 

This lack of overlap results partly from different dust de¬ 
tection limits between the two samples; the closer distance 
to HeViCS ETG and deeper observations result in much 


lower detectable dust levels in the HeViCS s ample. From the 
H-AT LAS Science Demonstration Phase, ([Rowlands et al.l 
I 2 OI 2 II stacked 233 optically selected ETG from GAMA and 
found that they have average dust masses of an order of 
magnitude less than their H-ATLAS detected ETG. How¬ 
ever, this does not explain why there are no HeViCS ETG 
occupying the same regions as the H-ATLAS sub-mm de¬ 
tected ETG. This cannot be explained by a detection limit, 
but may be due to the larger area surveyed by H-ATLAS. 
Some of the more extreme cases in H-ATLAS could be ex¬ 
plained as unusual ETG, but given that all H-ATLAS ETG 
occupy the top region of Fig. [SI this effect is likely explained 
by the difference in environment. Alternatively, this could 
be caused by the inability of dust to survive in dense envi¬ 
ronments due to galaxy-cluster interactions. Hydrodynam- 
ical or gravitational interactions that are likely to occur in 
dense regions may shorten the lifetime of dust, as may hot 
gas in the hostile intracluster medium. 

There are only three HeViCS ETG which demonstrate 
normalised dust masses on a level with the H-ATLAS ETG 
(VCC 327, 450 and 571). These all have GOLDMine clas¬ 
sifications of SO and are found in high density regions of 
Virgo (Egai ~100-200 gals Mpc“^). It may be possible that 
these galaxi es have been rece ntly accreted into the Virgo 
cluster ('e.g. [Kraft et al.ll201llf . and have not yet been sub¬ 
jected to the effects of dust stripping and destruction in the 
intra-cluster medium. 

We run a simple test to check whether this may be 
a possibility. Based on Virgo infall velocities provided by 
[Mould et al.[ ([ 2 OOOII an d assu ming a Virgo cluster radius 
of 2.2 Mpc ([Mei et al.[ [2007[l . we calculate typical cross¬ 
ing times for these three Virgo ETG of >^0.7, 0.8 and 2.2 
Gyr. If we assume ram-pressure stripping is responsible 
for the majority of dust loss in Virgo ETG, with typica l 
removal timescales of a few x 10® yr ([Takeda et al.[ [l984[ : 
[Murakami &: Babul[[l99^ . then it may be possible to relate 
these high relative dust levels with recent galactic infall into 
the Virgo cluster. These relatively short crossing times illus¬ 
trate that the current local density may not be the same as 
the time averaged local density experienced by a galaxy in 
the Virgo cluster. This is in contrast to what happens in the 
field. This may contribute to the different behaviours with 
environment shown in Fig. [5| 

Dust masses detected in these galaxies can be used to 
estimate total (both atomic and m olecular) gas mass: a typ¬ 
ical gas-to-dust ratio of 100 (e.g. [Parkin et aI.[[2(H^ gives 
a range of ~1x10^-8x10®Mq for the H-ATLAS sample and 
~lO®-5xlO®M0 in HeViGS. ATLAS®^ estimate molecular 
gas mass for some of the Virgo galaxies in this sample, find- 
ing an upper limit of 10 ® ®®Mq for these particular galaxies 
([Tbung et al.[ [201l[ ~l , which is consistent with our estima¬ 
tions for total gas mass. Additional results from ATLAS®^ 
find a strong HI detection rate dependence on surface den¬ 
sity whereby HI i n ETG is prefere ntially detected outside 
the Virgo Cluster ([Serra et al.[[2012f : again these results are 
qualitatively consistent with our findings for the two sam¬ 
ples whereby H-ATLAS ETG demonstrate a factor of ten 
dust and hence gas mass higher than HeViCS ETG in the 
dense regions of Virgo. However, while we note that dust 
and gas masses are both low in these Virgo ETG, S13 found 
that dust and H-I detections in Virgo ETG showed very 
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Figure 13. Distributions of Tform (left) and Tiastb (right) for H-ATLAS (red) and HeViCS (blue) populations. Numbers in the plots 
represent KS probabilities of parameter distributions coming from the same parent sample. Dashed lines represent mean values in the 
samples’ respective colours. Error bars represent the mean Icr errors for both samples. For Ti^stb (right) galaxies with NUV-r>5.0 are 
indicated with green histograms. 


small overlap, counter to the assumption that dust follows 
gas mass. 


Similar studies run on samples of LTG in the Virgo 
Cluster have demonstrated appreciably low levels of HI 
gas compared to LTG in sparser environments. Addi¬ 
tionally, lower star formation activity has been identi- 
hed in these spiral galaxies , and possibly lower dust lev¬ 
els teoselli fc Gavaz^ l2006l and references therein). Mod- 
els indicate that ram pressure stripping, gas compression 


els indicate that ram pressure stripping, g 
llBvrd fc Valtone^ll990l : lTonnesen fc Brvan 


20091) and star- 


__ 

vation due to the cluster potential ( Baloeh et al.ll20(^ are 
possible causes of these decreased levels of gas and dust in 
LTG. Theoretically, ETG in the same environment would 
also be subjected to these same physical mechanisms, re¬ 
sulting in the lower levels of gas and dust currently being 
observed. 


_ The Herschel Reference Survey (HRS, iBoselli et al.l 

l 2010 bl ) sampled a wider range of galaxy environments than 
just the Virgo cluster. Although it includes very few lumi¬ 
nous ETG other than those in the Virgo cluster, it is still 
useful to consider where their ETG reside in terms of param¬ 
eter space, a nd how this compare s to the two samples inves¬ 
tigated here. ISmith et al.l (120121 ) End 31 ETG in the HRS 
parent sample with 250^m detections: these all have stellar 
masses at > 1 O^°M 0 and modBB fits to the sample reveal 


a dust mass range of 1 O®'°“^'^M 0 and dust temperatures 
of 16-32K. Most notably, however, while they find a similar 
trend for dust-to-stellar mass ratio with stellar mass, their 
elliptical galaxies are found to present the lowest normalised 
dust masses. This is not what is seen here, particularly for 
the H-ATLAS/GAMA sample. It should be noted that the 
HRS sample only contains seven sub-mm detected elliptical 
galaxies, and therefore this result may be due to poor statis¬ 
tics. The majority (~ 68 %) of the HRS ETG reside within 
the Virgo cluster, which explains the similar dust mass ra nge 
to that of the HeViCS survey; in fact lSmith et al.l ll2012l) ex¬ 
plicitly state that there is overlap between their ETG and 
those of S13. Therefore we choose not to include a further 
study with HRS ETG. 

In summary we find tentative evidence that specific 
dust mass depends broadly on environment, however, more 
galaxy samples at intermediate environment density are 
needed to confirm such a trend. 


7 CONCLUSIONS 

This work has compared H-ATLAS sub-mm detected ETG 
to HeViCS (Virgo Cluster) sub-mm detected ETG. This was 
a strongly motivated study, as multiple Herschel works have 
revealed different levels of dust in different samples of ETG 
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dSkibba et al.|[20ll] : ISmith et al.ll20lll : [Rowlands et al.l|2012l : 

S13; A13). It has been unclear thus far whether these dif¬ 
ferences are simply due to different sample statistics and/or 
selection effects, or whether they are real differences which 
are a result of the different samples observing different types 
of ETG. 

Two samples were selected for this study: the A13 H- 
ATLAS sample of 220 ETG and 33 ETG from the HeViCS 
S13 sample, both with Mr < —17.4mag. With the aid of con¬ 
sistent calculations for nearest neighbour density, and MAG- 
PHYS panchromatic SED fitting to the multi-wavelength 
data, we were able to objectively quantify the true differ¬ 
ences in the properties of these ETG. The results of this 
study are summarised here below. 


(i) Nearest neighbour surface densities revealed true dif¬ 
ferences in the type of environment in which these ETG re¬ 
side. H-ATLAS ETG are in isolated environments, spanning 
<0.1-10 galaxies Mpc“^, whereas HeViCS ETG are domi¬ 
nated by the cluster environment ('^25-500 galaxies Mpc“^). 
These results are also true for undetected ETG in each sam¬ 
ple, with only a trace overlap in density between samples 
observed at ~20-100 galaxies Mpc“^. We find that sub-mm 
detected ETG in H-ATLAS reside in sparser environments 
than undetected ETG. 

(ii) ModBB fits from A13 and S13 reveal different ranges 
of dust-to-stellar mass ratio, with H-ATLAS ETG demon¬ 
strating higher Md/M« at fixed stellar mass. We prove that 
this is not a selection effect by carrying out a KS test and 
by Monte Carlo simulations of Virgo ETG at GAMA ETG 
distance. Both these tests confirm that the samples have 
< 0.1% probability of having been drawn from the same 
specific dust mass distribution. 

(iii) MAGPHYS results indicate that it is sometimes diffi¬ 
cult to accurately constrain the cold dust temperature, but 
similar results for dust mass were obtained using MAG¬ 
PHYS and ModBB fits, in spite of these uncertainties. 
ModBB fits appear to give higher dust mass of a galaxy 
in some cases - this may be because the smaller dust grains 
which emit at higher temperatures are not accounted for in 
the ModBB fits. 

(iv) MAGPHYS stellar masses incl uding UKIDSS fluxes 
are lower than those produced by Zibetti et al. ( 2009l l 
(HeViCS) and the GAMA team (iTavlor et al. 201] ) (H- 
ATLAS). Both sets of stellar masses indicate that H-ATLAS 
ETG are more massive on average than HeViCS ETG. 

(v) Correlations are found between dust mass and stellar 
mass for both H-ATLAS (rp=0.42) and HeViCS (rp=0.58) 
ETG. Additionally anticorrelations are found between dust- 
to-stellar mass ratio and stellar mass, although the trend 
is shifted upwards (to higher normalised dust mass) for H- 
ATLAS. Most of this anticorrelation is due to the dust de¬ 
tection limits. However it remains to be understood why 
there is a lack of massive ETG in Virgo with high dust-to- 
stellar mass ratios. Investigating dust-to-stellar mass ratio as 
a function of nearest neighbour density reveals another cor¬ 
relation between the two properties, where both H-ATLAS 
and HeViCS ETG sit on the same trendline. This is an indi¬ 
cator that levels of dust mass in ETG are affected by their 
environments. 

(vi) Examinations of sSFR reveal that dust mass is indica¬ 
tive of ongoing star formation in many of these galaxies. 


but is not directly related, as evidenced by different trends 
in specific dust mass and sSFR plots with stellar mass. It 
appears that there is very little (if any) ongoing star forma¬ 
tion in the HeViCS ETG, but quite the opposite is true for a 
large proportion of the H-ATLAS sample, with the highest 
sSFRs on par with that of spiral galaxies. 

(vii) The massive ETG in H-ATLAS have similar sSFRs to 
the HeViCS ETG. 
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Figure Al. Example galaxy VCC 1535 in all four WISE bands. 
In each frame the blue ellipse is the parent aperture for that band 
from which concentric ellipses were derived. In the W1 frame (top 
left), the masked foreground stars are indicated in red circles. 
These pixels were masked in all bands. 


APPENDIX A: WISE PHOTOMETRY 

The Wide-field Infrared Survey Explorer (WISE) has a beam 
size of 12.0" (in the W4 band) as its coarsest resolution. 
This means most of the galaxies in the Virgo cluster are 
extended objects in all WISE bands. Magnitudes from the 
WISE Science Archive 0 are based on aperture photometry 
using an elliptical aperture for each galaxy. Eoreground stars 
are not removed, hence diluting the flux measurements. We 
therefore opted to compute asymptotic fluxes for all galaxies 
in our sample, taking special care of the contamination by 
bright foreground stars. 

The images for each galaxy were retrieved from the 
WISE Science Archive. We choose sufficiently large cut-outs 
from the All-Sky Atlas in order to cover both the galaxy and 
its close surroundings. 

The local background of the galaxy is estimated from 10 
boxes of 50x 50 pixels. Three values are computed from these 
boxes: 1) a single background value. A sigma clipped mean 
was derived from all pixels using a 5cr rejection threshold and 
iterations until convergence. The background value was then 
subtracted from the image. 2) A pixel-to-pixel background 
error. This is the mean of the standard deviations on the 
pixel values in each of the boxes. 3) A large scale background 
variation error. This is the standard deviation of the mean 
values of each of the boxes. 

We identify the brightest stars and mask their corre¬ 
sponding pixels so they will be ignored in the further analy¬ 
sis. An elliptical aperture is then determined, following the 
galaxy’s apparent shape on the sky. The angle and ellipticity 
of this parent aperture will serve as a base for constructing 


Science Archive: http://irsa.ipac.caltech.edu 




Figure A2. Top: Growth curve for the W1 band of VCC 1535, 
showing the log of the cumulative counts for each elliptical shell. 
The black vertical line indicates the major axis of the parent aper¬ 
ture. Bottom: Cumulative counts as a function of the gradient 
of the growth curve. The red line is a linear fit (displayed here 
on log-linear scale) to the blue points in that interval. The green 
line is the extrapolation of this linear fit to all of the points. 


the growth curve. To construct the growth curve for each 
galaxy and each band, concentric elliptical annuli are cre¬ 
ated based on the parent aperture. Starting from the center 
of the galaxy, each next aperture is a fraction An in ma¬ 
jor axis larger. We set Aa = a/20, where a is the major 
axis of the parent aperture. The total flux inside the ellipti¬ 
cal shell is estimated by multiplying the median value of all 
pixels inside the shell with the total number of pixels. In the 
number of pixels we also include the amount of masked pix¬ 
els. This extra step filters out any remaining contribution of 
foreground stars, while approximating the flux of the galaxy 
in that shell. As soon as the calculated flux in a shell falls 
below the background level, the iteration is stopped. 

The growth curve is now plotted as the cumulative 
counts ve rsus the distance to the c enter of the galaxy. 
Following iMufioz-Mateos et al.l ll2009l) , the gradient of the 
growth curve around the edge of the galaxy is calculated. 
Fig. IA2I (lower panel) shows the cumulative counts as a 
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function of this gradient. A linear fit to these points is per¬ 
formed and the intercept with the y-axis is calculated as 
the asymptotic flux of the galaxy in that particular band. 
As the gradient changes rapidly and is non-linear inside the 
parent annulus of the galaxy, these points are not consid¬ 
ered in the fit. The outer points of the growth curve are not 
considered in the fit either, as they are usually dominated 
by background variations or unmasked foreground stars. 

Four sources of uncertainty on the calculated fluxes are 
considered, a) The Poisson noise, determined as the square 
root of the asymptotic flux. A multiplicative factor of 2 for 
W1,W2 and W3 and 4 for W4 must be applied to incorpo¬ 
rate the correlated pixel noise, according to the Explanatory 
Supplement to the WISE Preliminary Data Release Prod¬ 
ucts 0, section II.3.i. h) The pixel-to-pixel background er¬ 
ror as described above, c) Large scale background variations, 
also described above, d) A calibration uncerta inty. Following 
the recommendations of ijarrett et al.l il201ll L we use 2.4%, 
2.8%, 4.5% and 5.7% as the calibration error in Wl, W2, 
W3 and W4, respectively. 

The obtained fluxes and their respective errors have to 
be converted to Jy to be of physical meaning. Section II.3.f 
of the WISE photometry manual PI. provides the following 
conversion factors: 1.9350 x 10“® Jy/DN , 2.7048 x 10“® 
Jy/DN, 2.9045 x 10"® Jy/DN and 5.2269 x 10"® Jy/DN for 
Wl, W2, W3 and W4, respectively. 

Ijarrett et al.l (l2013lf advises three corrections to the 
WISE flux of extended sources. One of them is a colour 
correction, which we do not apply here as our fitting rou¬ 
tine takes the filter response into account. The second cor¬ 
rection stems from a calibration d iscrepancy between blue 
stars and red galaxi es, described bv lWright et al.l (l2010ll and 
Ijarrett et al.l 1 I 2 OIII L They advise a multiplicative factor of 
0.92 for the W4 flux of all spiral and disk galaxies. As we 
are dealing with early type galaxies in our sample, we have 
no need of such a correction. This leaves us with the third 
correction, which is an aperture correction due to the fact 
that the absolution calibration for WISE was done using 
PSF profile fitting. We apply a correction of 0.03 mag, 0.04 
mag, 0.03 mag and —0.03 mag for Wl, W2, W3 and W4, 
respectively. 

We verify our measurements by comparing them with 
the archival and literature WISE fluxes, determined by flux 
measurements of elliptical apertures (see Fig. IA3I) . The 
slopes of the data points in each of the panels are roughly 
parallel to the 1 : 1 relation (black line). There is, how¬ 
ever, an offset visible in all bands when comparing to the 
archival fluxes (blue points). Interestingly, we find system¬ 
atically higher fluxes in the Wl and W4 band, while our 
mea surements yield low er fluxes in the W2 and W3 bands. 

ICiesla et al.l ll2014l ') measured the W3 and W4 fluxes of a 
small s ubset of our sample a s part of the Herschel Reference 
Survey llBoselli et al.l2010al L The corresponding galaxies are 
of the largest in our sample and include the 4 radio galax¬ 
ies. Although the measurements were also done using aper¬ 
ture photometry, much more care was taken in the choice of 
the apertures than the automated WISE pipeline. Further¬ 
more, the contribution of foreground stars is less significant 





Figure A3. Comparison of the integrated WISE fluxes as derived 
by our method (asymtotic fluxes) and the archival fluxes (blue 
points) measured from elliptical annuli. Red points are indepen¬ 
dent p hotometric measurements in W3 and W4 from ICiesla et aO 
1120141 1 for a subset of our sample. The black line is the 1 : 1 rela¬ 
tion. 


in these two b ands. Our agreement in both W3 and W4 with 
ICiesla et al.l ll2014l ) boosts confidence in our method and 
recommends independent flux measurements over archival 
WISE fluxes when it comes to extended sources. Our re¬ 
sultant fluxes for these WISE bands are displayed fully in 
Table lAll 


http://wise2.ipac.caltech.edu/docs/release/prelim/expsup/wise_prelrel_toc.html 
See footnote 1101 
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Table Al. WISE photometry for the HeViCS sample. All measurements are in mjy. 


vcc 

-Pwi 

F'w2 

Fw3 

Fw4 

94 

42.8 ± 1.1 

23.54 ± 0.78 

22.9 ± 1.1 

- 

220 

100.3 ± 2.6 

51.6 ± 1.6 

26.9 ± 1.3 

13.4 ± 1.8 

270 

8.12 ± 0.27 

3.54 ± 0.18 

8.89 ± 0.49 

12 ± 1.8 

278 

6.58 ± 0.23 

3.73 ± 0.18 

- 

- 

312 

51.8 ± 1.4 

26 ± 0.86 

13.02 ± 0.68 

- 

327 

4.18 ± 0.17 

2.52 ± 0.14 

3.77 ± 0.24 

4.9 ± 1.1 

345 

515 ± 13 

266.6 ± 7.9 

110.2 ± 5.2 

74.1 ± 5.7 

408 

233.5 ± 5.9 

117.6 ± 3.5 

85.8 ± 4.1 

56.8 ± 4.7 

411 

18.61 ± 0.54 

10.4 ± 0.4 

- 

- 

450 

5.3 ± 0.2 

2.16 ± 0.13 

8.72 ± 0.47 

10.8 ± 1.6 

462 

64 ± 1.7 

42.4 ± 1.3 

29.1 ± 1.4 

16 ± 2 

482 

13.2 ± 0.4 

8.72 ± 0.34 

3.67 ± 0.25 

- 

571 

5.4 ± 0.2 

2.97 ± 0.16 

3.43 ± 0.23 

5.7 ± 1.2 

672 

17.54 ± 0.51 

10.6 ± 0.4 

7.78 ± 0.43 

10.6 ± 1.6 

685 

245 ± 6.2 

134 ± 4 

67.7 ± 3.2 

34.7 ± 3.3 

758 

74.3 ± 1.9 

39 ± 1.2 

- 

- 

763 

1417 ± 35 

775 ± 23 

316 ± 15 

92.4 ± 6.8 

781 

7.29 ± 0.25 

3.88 ± 0.19 

2.9 ± 0.2 

- 

881 

1902 ± 47 

1026 ± 30 

298 ± 14 

79 ± 6 

951 

13.3 ± 0.4 

6.06 ± 0.26 

3.69 ± 0.24 

- 

1003 

697 ± 17 

370 ± 11 

222 ± 10 

103.9 ± 7.5 

1030 

311.8 ± 7.9 

169.5 ± 5.1 

114.1 ± 5.4 

123.3 ± 8.6 

1154 

553 ± 14 

274.5 ± 8.2 

215 ± 10 

130 ± 9 

1226 

3203 ± 80 

1307 ± 38 

587 ± 28 

150 ± 10 

1250 

59 ± 1.5 

33.1 ± 1.1 

43.6 ± 2.1 

35.4 ± 3.4 

1253 

430 ± 11 

225.7 ± 6.7 

79.4 ± 3.8 

38.7 ± 3.6 

1316 

2275 ± 57 

1067 ± 31 

464 ± 22 

208 ± 13 

1535 

854 ± 21 

437 ± 13 

407 ± 19 

317 ± 20 

1614 

10.5 ± 0.33 

6.82 ± 0.28 

3.65 ± 0.24 

1.82 ± 0.63 

1619 

112 ± 2.9 

57.9 ± 1.8 

23.3 ± 1.2 

9 ± 1.5 

1632 

922 ± 23 

451 ± 13 

160.5 ± 7.6 

80.2 ± 6.2 

486 

17 ± 0.5 

10.19 ± 0.39 

2.9 ± 0.2 

5.6 ± 1.1 

1327 

56.5 ± 1.5 

78.6 ± 2.4 

15.53 ± 0.79 

9.5 ± 1.5 






